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ABSTRACT 

We present the discovery of a high redshift, X-ray selected AGN by the Calan- 
Yale Deep Extragalactic Research (CYDER) survey: CXOCY J033716.7-050153, 
located at 2; = 4.61, the second high redshift AGN discovered by this survey. 
Here, we present its optical, near-IR and X-ray properties and compare it with 
other optical and X-ray selected high redshift AGN. The optical luminosity of this 
object is significantly lower than most optically selected high redshift quasars. 
It also has a lower rest frame UV to X-ray emission ratio than most known 
quasars at this redshift. This mild deviation can be explained either by dust 
obscuring the UV radiation of a normal radio quiet AGN emitting at 10% of 
its Eddington luminosity or because this is intrinsically a low luminosity radio 
loud AGN, with a super-massive Black Hole of ~ IO^Mq emitting at 1% of its 
Eddington luminosity. Deep radio observations can discriminate between these 
two hypotheses. 

Subject headings: galaxies: active — quasars: individual (CXOCY J033716.7- 
050153) — X-rays: galaxies 
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1. Introduction 

Due to its large intrinsic luminosity, an AGN (Active Galactic Nucleus) can be observed 
up to very high redshifts, and AGNs can therefore be used to probe the early epochs of the 
universe when the formation of large structures began. Recent observations have found a 
considerable number of quasars at high redshift (Anderson et al. 2001; Schneider et al. 2002), 
with some of them at redshift > 6 (Fan et al. 2003), when the universe was about 1 Gyr old. 

While AGNs were first noticed to be a class of objects worthy of detailed follow-up in 
other bands on the basis of their radio emission, e.g. (Schmidt 1963), more recent work has 
focused on surveys at shorter wavelengths. Due to their low spatial density, most AGNs at 
high redshift (defined here to be z > 4) were discovered by large area, shallow optical surveys. 
An example of this is the Sloan Digital Sky Survey (SDSS; York et al. 2000), which in its 
first data release (Schneider et al. 2003) presented 238 high redshift quasars. However, these 
quasars were selected based on their optical properties, and optical surveys are subject to 
bias against the discovery of quasars which are either intrinsically dim at these wavelengths 
or are obscured by dust like the type 2 AGN that are believed to comprise most of the AGN 
background (Hasinger 2002). 

The other main AGN selection technique is X-ray detection, as X-ray emission appears 
to be a universal characteristic of AGNs at all observed redshifts (Kaspi et al. 2000). How- 
ever, given the relative faintness of the majority of the sources. X-ray surveys need to rely 
on optical spectroscopy for most of the identifications. With new X-ray missions like the 
Chandra X-ray Observatory and its superb angular resolution, the identification of X-ray 
sources at other wavelengths has become much easier, allowing for successful multiwave- 
length follow-up of X-ray sources. The nature of the X-ray emission also ensures that the 
optical selection bias due to obscuration is strongly reduced (although the X-ray selection 
may introduce other biases). This is especially true at high redshift, where the observed 
X-ray photons were emitted at higher energies and can penetrate even considerable amounts 
of obscuring material. 

According to a recent review by Brandt et al. (2003) only seven of the high redshift 
AGNs known were discovered by their X-ray emission. Three of them were found with 
ROSAT (Henry et al. 1994; Zickgraf et al. 1997; Schneider et al. 1998), while the other four 
were detected by Chandra (Silverman et al. 2002; Barger et al. 2002; Castander et al. 2003a). 
One of them was detected in the radio and classified as radio loud (Zickgraf et al. 1997). 

In this paper, we present the discovery of the eighth high redshift X-ray selected AGN, 
and the second discovered by the Calan Yale Deep Extragalactic Research (CYDER) survey, 
at a redshift z = 4.61. The CYDER survey is a collaborative effort between Universidad de 
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Chile and Yale University to study faint stellar and extragalactic populations in detail. One 

of the key aims is the characterization of the population of faint X-ray sources. In order to 
do this, some of the fields of the survey were selected to overlap with moderately deep (~ 50 
ks) Chandra pointings. Multiwavelength follow-up in these fields includes optical imaging 
using the CTIO 4m telescope, near IR images obtained using the DuPont telescope at Las 
Campanas Observatory (LCO) and optical spectroscopy from VLT at Ccrro Paranal and 
Magellan at LCO. A more complete review of the optical and spectroscopic properties of the 
X-ray sources in the first two fields can be found in Castander et al. (2003b). 

In §2 wc describe the X-ray, optical and ncar-IR observations that led to the discovery 
of the AGN presented here. In §3 we discuss its observed properties and compare it with 
other high redshift quasars, both optical and X-ray selected. We present our conclusions 
in §4. Throughout this paper we assume Hq — TOhm km s Mpc~^, = 0.3 and 
= 0.7. These cosmological parameters are consistent with the recent findings by Spergel 
et al. (2003). Wc define the photon index F as the exponent giving a photon fiux density 
X-ray spectrum dN/dE oc in photons cm~^ s~^ keV~^. 

2. Observations and Data Analysis 
2.1. X-Ray Data 

As part of the CYDER Survey, five fields observed by Chandra and currently available 
in the archive were downloaded and analyzed using standard techniques with the CIAO 
package. In this paper wc discuss the results for the source CXOCY J033716. 7-050153 found 
in the SBS 0335-05 field. This field was observed by Chandra on September 7th, 2000 for 
60.51 ks (PI: Thuan) and retrieved from the archive by our group. A detailed analysis of 
the whole sample of X-ray sources detected in these fields will be presented in Treister et al. 
(2003). 

Reduction of the data included the removal of bad columns and pixels using the guide- 
lines specified on the "ACIS Recipes: Clean the Data" web page and the removal of flar- 
ing pixels using the FLAGFLARE routine. We used the full set of standard event grades 
(0,2,3,4,6) and created two images, one from 0.5 to 2.0 keV and one from 2.0 to 8.0 keV. 
Then, we used the WAVDETECT routine from the CIAO package to identify the point 
sources within these images, checking wavelet scales 1,2,4,8 and 16. 

CXOCY J033716.7-050153 was detected in the 12 ACIS CCD. We extract an X-ray 
spectrum using a minimum of 5 counts per channel from 0.5 to 3.5 keV (i.e. the 20 lowest 
energy photons), using the PSEXTRACT script. We then fltted the spectrum within XSPEC 



-4- 



11.0, with a power law model and absorption constrained to the Galactic value of 4.8 x 10^° 
from the Nh tool. This is relatively unconstraining, yielding a spectral index of F = 2.0±0.6, 
and observed fluxes of 1.5 x 10^^"'' ergs/sec/cm^ and 3.2 x 10^^"'' ergs/sec/cm^, from 0.5-2 
and 0.5-8 keV, respectively. The fit to the X-ray spectrum did not depend on whether the 
or the C-statistic was used, so it is unlikely that large errors were induced by fitting 
with an insufficient number of photons per bin. The unabsorbed fluxes are approximately 
1.8 X 10~^^ and 3.5 x lO"-*^^ ergs/sec/cm^, respectively. The errors should be about 30% due 
to the statistical error in the count rate and the uncertainties in the best-fitting spectral 
models. The cumulative histogram of detected X-ray photons is presented on Figure 1. 

2.2. Optical Imaging Data 

The field of the AGN was observed with the CTIO 4 meter Blanco telescope, using the 
MOSAIC-II camera, that provides a field of view of 36'x36'. Images were taken using the 
By,R and / filters on the nights of October 8th and 10th, 2002, centered on the Chandra 
archival image central coordinates. Exposure times were 2 hours under 1.3" seeing conditions 
in the B filter, 50 minutes under 1.1" seeing in V, 50 minutes under 1.1" seeing in R and 20 
minutes under 1.0" seeing in I. 

A second image in the / band was taken with the UT4 VLT telescope at Cerro Paranal, 
Chile using the F0RS2 instrument. This image was taken on the night of September 19 
2003 in service mode. Total exposure time was 14 minutes under 0.7" seeing conditions. 

Images were reduced using standard techniques with the IRAF/MSCRED package. 
Identification of the optical counterpart of the X-ray source was straightforward given the 
superb spatial resolution of the Chandra Observatory. Indeed, only one optical source can 
be found in a 3" circle around the X-ray emission centroid. The offset between the X-ray 
and optical positions was < 0.1" in RA and 0.5" in Dec. In order to calculate the received 
fiux in each band, we performed aperture photometry using an aperture of 1.4xFWHM of 
each image, centered on the centroid of the R band source, the band in which the signal to 
noise is highest. The observed magnitudes in each band can be found in Table 1. Figure 2 
shows the optical finding chart, based on the R band image. 

2.3. Optical Spectroscopy 

Multislit spectroscopy of X-ray sources detected by Chandra in the field of the AGN was 
obtained at Las Campanas Observatory with the Magellan I (Baade) telescope, using the 
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LDSS-2 instrument. The observations were taken on the night of October 4th, 2002. The field 
of the high redshift AGN CXOCY J033716. 7-050 153 was observed for 2 hours in 0.75" seeing 
conditions. The Med/Bluc grism was used, giving a dispersion of 5.3A per pixel at a central 
wavelength of 5500A. The obtained spectrum was reduced using standard IRAF tasks called 
by a customized version of the BOGUS code^ .We calibrated the wavelength of the spectrum 
using the He-Ar comparison lamp and the night sky lines. Rough flux calibration of the 
spectrum was performed using the spectrum of the LTT9239 spectrophotometric standard. 

Figure 3 shows the final, reduced spectrum together with the mean spectrum of quasars 
at z > A observed by the SDSS. Lya and NV arc clearly distinguishable in the spectrum. 
The measured redshift using these emission lines is z = 4.61 ± 0.01. 



3. Discussion 

With an absolute magnitude of Mb = —21.14 (Vega), calculated extrapolating the 
I band magnitude using the SDSS QSO composite spectrum (Vanden Berk et al. 2001), 
CXOGY J033716. 7-050153 is a faint AGN. Indeed, there arc only two AGNs known at high 
redshift that are fainter in the optical, VIA J1236-F6213 at z = 4.42 (Brandt et al. 2001) and 
CXOHDFN J123719.0+621025 at z ^ 4.13 (Barger et al. 2002), both located in the HDF- 
N/CDF-N. In X-rays, this quasar shows a luminosity in the [0.5-2.0] keV band of 3.77 x 10^^ 
ergs s~^. This is a large X-ray luminosity for an AGN that is faint in the optical bands. In 
order to quantify this statement, we can calculate the effective optical to X-ray power-law 
spectral slope which is given by 

_ log[/.(2 keV)//,(2500 A]] 
"^""^ \og[u{2 keV)/i/(2500 A)] ^ ^ 

where is the fiux density per unit of frequency and u is the frequency of the given wave- 
length or energy. For this AGN, we measure a value of = — 1.16 ± 0.16. The average 
value for the known sample of 2; > 4 AGNs detected in X-rays^ is —1.61 ± 0.24 from a total 
of 71 sources, so this new AGN is a ~ 2(7 deviation. With the data available now we cannot 
calculate the contribution from the host galaxy to the optical luminosity, but it can only be 
significant if its luminosity is > 3L*. Furthermore, a significant contribution from the host 
galaxy will only make the optical luminosity of the AGN to be smaller, therefore increasing 

''Available at littp://zwol£k;inder.jpl.nasa.gov/'~stern/homepage/bogus.html 

*^Tliis sample can be found in the Web site http://www.astro.psu.edu/user/mel/papers/liigliz-xray- 
detected.dat maintained by Niel Brandt and Christian Vignali 
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the value of aox and making the deviation from a typical AGN at this redshift to be even 
larger. 

A plot of function of redshift for all the high redshift AGNs detected in X-rays 

can be found in Figure 4 and the Fx / Fuv relation for these objects is presented in Figure 
5. In Figure 4 we can see that the X-ray selected sources have systematically higher values 
of aox than optically selected AGNs. This correlation is expected from selection effects; 
however, it should be noted that the optical luminosities of all the X-ray selected sources are 
fainter than that of the least luminous optically-selected quasar (with the exception of the 
one X-ray selected blazar), see Figure 5. 

Recents results presented by Vignali et al. (2003b) show that there is a good correlation 
between the UV flux characterized by the flux density at 2500A and the soft X-ray flux 
measured at 2 keV. This correlation can be expressed as Fx oc F^'^ . A similar result was 
reported in (Castander et al. 2003a), in this case Fx oc F^j-y'. This relation implies that aox 
should also be correlated with the UV flux (or luminosity). Figure 5 shows that optically 
selected AGNs are more luminous, both in the optical and the X-ray, than the X-ray selected 
ones. The agx value of optically selected AGN are therefore expected to be lower as can be 
seen in Figure 4. 

In Figure 5, we can also see that without correcting for the possible optical obscuration 
CXOCY J033716.7-050153 is a significant deviant from the Fx - Fuv relation. There are 
two possible explanations for this deviation: either there is significant obscuration of the 
optical light by dust or this a radio loud AGN with aox enhanced by a second component, 
probably coronal, contributing to the X-ray emission. Strong optical variability can be ruled 
out using the two optical observations separated by ~ 1 year in the observed frame, in which 
the / band magnitude remained almost unchanged given the error bars. 

If this AGN is relatively unobscured, then its anomalous value of aox-, as well is its 
relative optical faintess compared with other high redshift AGN may be explained if the 
source is at the bright end of the low luminosity AGN (LLAGN) category. Such sources are 
typically seen to have X-ray luminosities below about 1% of the Eddington limit, e.g. Ho 
(2002); Maccarone et al. (2003). The mass of the central black hole would then have to be at 
least about 3.4 x 10^ M©. Applying the black hole fundamental plane relation from Merloni 
et al. (2003), 

log Lr = (0.601°:}}) log Lx + (0.781°:^^) log M + 7.331^2^ (2) 

where Lr is the radio luminosity at 5 GHz, Lx is the 2-10 keV X-ray luminosity and M is 
the mass of the black hole in units of Mq, we find that the radio power is expected to be at 
least fsGHz — 31.5/xJy. Now, if we assume instead that this is a typical AGN radiating at 
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about 10% of its Eddington luminosity then using the Merloni relation we get a radio flux of 
fbGHz = 5.2/iJy. Moreover it seems that the radio emission is quenched in AGN accreting at 
about 10% of the Eddington rate (Maccarone et al. 2003), so an even lower radio flux would 
be expected. 

Following the standard definition of radio to optical flux ratio Rro given by Kellermann 

et al. (1989), where optical flux is defined as the flux at 4400A in the rest frame and radio 
flux is calculated at a rest frame frequency of 5 GHz (6 cm), radio loud sources have Rro 
values in the range 10-1000, while for radio quiet AGNs 0.1< Rro <1- With the assumption 
that the Eddington ratio is 1% for this AGN we get that Rro = 15.06, making it a radio 
loud AGN. However, given the low optical luminosity, we cannot neglect the contribution of 
star light from the host galaxy, in which case the value of Rro computed is simply the lower 
limit of the radio to optical flux ratio for the AGN emission alone. Assuming an Eddington 
ratio of 10% then Rro — 2.49, which makes it a radio quiet AGN. 

According to an alternate definition of radio loudness often used for high redshift quasars 
(Stern et al. 2000), to be considered radio loud an AGN has to have a minimum 1.4 GHz 
specific luminosity of L-^ 4 qjj^ — 1-61 x 10^^ ergs s~^ Hz~^ (Schneider et al. 1992), which at 
z — 4.61 corresponds to a specific flux of 4 — 0.9 mJy, very close to the detection limit 
of the VIA FIRST survey. Unfortunately, the field of CXOCY J033716. 7-050153 was not 
observed by the FIRST survey. Using the NVSS survey (Condon et al. 1998) catalog, there 
were no detectable sources within a 1 arcminute radius. Considering that the completeness 
limit is 2.5 mJy, we can use this as an upper limit to the 1.4 GHz continuum emission. Then, 
this AGN could be at most marginally classified as radio loud. 

In order to determine the amount of obscuration affecting the optical spectrum of this 
AGN we can compare its optical colors with a typical AGN spectrum at this redshift. The 
optical colors of this AGN arc: V - R = 1.75 ± 0.40 and V - I = 2.1 ± 0.39. Using the 
composite SDSS QSO spectrum (Vanden Berk et al. 2001) and redshifting it to 2; = 4.61 
we obtain the following colors: V — R — 1.02 and V — I — 1.71 accounting for the IGM 
absorption using the description given in Madau (1995). If instead we use an average of all 
SDSS QSOs ai z > A the results are very similar, finding differences < 10%. The fact that 
CXOCY J033716. 7-050153 is shghtly redder than the average QSO at that redshift suggests 
that the optical spectrum is subject to some obscuration in the line of sight. In order to 
investigate this hypothesis, we added dust obscuration to the average spectrum following the 
prescription given in Cardelli et al. (1989). Given that the observed R band flux is affected 
by the Lya emission line, whose strength is highly variable from AGN to AGN, we will base 
our reddening estimate only on the V — I color. If order to obtain a, V — I color of 2.1 
for the redshifted composite SDSS quasar, similar to the observed value we need to add an 
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extinction contribution of Ay — 0.41 ± 0.4 in the rest frame. Therefore, we will adopt this 
value to correct the observed optical fluxes for intrinsic dust extinction. It is also worth 
noting that the accretion disk may become cooler at lower luminosities, which would also 
yield redder colors. 

Another typical way to calculate the amount of obscuration in the line of sight is using 
the X-ray spectrum. In this case, however, we do not have enough counts in the X-ray 
spectrum to calculate absorption directly. Also, since the observed soft X-ray band [0.5- 
2.0] keV translates into [2.8-11.22] keV in the rest frame, the observed X-ray spectrum is 
insensitive to moderate amounts of neutral hydrogen absorption {Nh < 10^^ cm~^). Using 
the standard dust to gas ratio, we can convert the optical extinction to a neutral hydrogen 
column density, given by Nh = 1-96 x lO^^^y cm~^ (Granato & Danese 1994). Then, we 
obtain in this case that Nh = 8 x 10^° cm~^, implying that the X-ray spectrum is not affected 
by absorption in the observed energy range. 

Now, correcting for the optical obscuration we compute an unobscured / magnitude of 
22.52±1.1. Using this value, we calculate the unobscured optical to X-ray power law slope 
— —1.35, which is now similar to the average value found for high redshift AGN. 

As we can see in Figure 5 and already presented in Castander et al. (2003a), fx oc f^y is 
a better fit to the data if the low UV luminosity X-ray selected AGNs are included. The values 
for CXOCY J033716. 7-050153 are consistent with this correlation if optical obscuration is 
taken into account. 



4. Conclusions 

We present the discovery of the CXOCY J033716. 7-050153, the second high redshift 
X-ray selected AGN discovered by the CYDER survey. This is a faint object in the optical. 
In fact there are only two AGNs fainter at these high redshifts. The X-ray fiux relative to 
its optical emission is unusual, having a larger value of than the typical AGN at high 
redshift. There are two effects that could bring this object into better accordance with the 
observed correlations for the general high redshift population: 

The first possibility is that this is a radio loud low luminosity AGN with a Black Hole 
of ~ 3 X 10^ Mq , emitting at ~ 1% of its Eddington Luminosity. In this case the accretion 
disk is cooler than a typical AGN at high redshift, making the optical colors look redder. 
Also, the high value of aox, even taking into account the correlation with optical luminosity, 
can be explained by the presence of second emission component, probably coronal, which 
contributes to the X-ray flux measurement. 
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A second possibility is that this is a radio quiet AGN emitting at an Eddington ratio 
of ~ 10%, in which the optical emisssion is obscured by a moderate amount of dust which 
would imply a neutral hydrogen column density of ~ 10^^ cm^^, causing the observed optical 
colors to be redder than the normal AGN. If this is the case, the observed value of is 
affected by the obscuration of the optical light, and after correcting for it, the corrected 
value agrees very well with the measurements for other high redshift AGNs, including the 
observed correlation between fx and fuy. 

Detection of a fiat radio spectrum in this source above about 10 yuJy would confirm the 
Low Luminosity AGN hypothesis, while a non-detection in the radio with a sensitivity limit 
of a few /xJy would strongly support the radio quiet, mildly obscured AGN hypothesis. 

We like to thank Meg Urry for useful discussions and the anonymous referee for his 
comments. The CYDER survey participants and especially FJC, ET and EG acknowledge 
support from Fundacion Andes. JM gratefully acknowledges support from the Chilean Cen- 
tro de Astrofisica FONDAP 15010003. EG acknowledges the support of an NSF Astronomy 
& Astrophysics Postdoctoral Fellowship under award AST 02-01667. 
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Fig. 1. — Cumulative histogram of detected X-ray photons. The sohd hne shows the distri- 
bution of the data, while the dashed line shows the fitted model with Nh — 4.8 x 10^° cm~^ 
and r = 2.0. Errors in the measurement are Poissonian. 
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Fig. 2.— R band image of CXOCY J033716.7-050153. This image is a cutout of 30"x30" of 
the original image, centered on the position of the optical counterpart of the X-ray source. 
North is up and east is to the left. This image was taken at the CTIO 4m telescope using the 
MOSAIC-II camera. A 2" circle centered on the optical counterpart position is also plotted. 
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Fig. 3. — Observed optical spectrum of CXOCY J033716. 7-050153 shifted into the rest frame 
using the best fit redshift of 4.61. Original dispersion is 5. 3A per pixel. The spectrum shown 
was boxcar smoothed using a 5 pixel box. Resolution of the rest frame smoothed spectrum 
is 5.6A . For comparison we also show the error- weighted average of the SDSS Early Data 
Release QSO spectra ai z > 4. See (Castander et al. 2003a) for details. 
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Fig. 4. — Plot of vs. redshift for the known sample of high redshift AGN detected 
in X-rays. Filled circles are ROSAT detections, empty circles are Chandra detections and 
the filled triangle is the only XMM-Newton detection. Stars are radio loud AGNs. Symbols 
enclosed by squares denote X-ray selection, those enclosed by triangles denote radio selection, 
while no enclosing sign denotes optical selection. The horizontal line marks the position of 
= —1.61, the average of the agx distribution. The large circle encloses the position of 
CXOCY J033716.7-050153. 
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Fig. 5. — /2]^ vs /25QQ ^ (^s defined in the text) for tlie known sample of liigh redsliift 
AGN detected in X-rays. Symbols arc the same as in Fig 4. Segmented lines show the 
positions of the ]x f^y^ and fx oc fjj'y' correlations reported by (Vignali et al. 2003b) 
and (Castander et al. 2003a) respectively. The long-dashed line shows the position of sources 
with fx oc fuv and aox = —1.61. The cross enclosed in a square marks the position of 
CXOCY J033716. 7-050153 after correcting for optical obscuration, while the arrow indicates 
the displacement caused by this correction. The uncertainty of the X-ray flux measurement 
of CXOCY J033716.7-050153 is ~ 30% and is not included for clarity. 



-16- 



Table 1. Properties of CXOCY J033716.7-050153 



Parameter 


Value 


RA (J2000) (X-ray) 


03''37™16.6s 


Dec (J2000) (X-ray) 


-05°01'54.3" 


RA (J2000) (optical) 


03''37"*16.6" 


Dec (J2000) (optical) 


-05°01'53.7" 


Galactic Nh^ 


4.8 X lO^o cm-= 


z 


4.61 ±0.01 




-1.16 ±0.16 


B mag^ 


> 26.5 


V mag*^ 


25.76±0.39 


R mag*^ 


24.01±0.08 


I mag (CTIO)^ 


23.95±0.23 


I mag (VLT)b 


23.66±0.07 


Mf 


-22.14±0.15 


^-Sl450(l+z) 


23.82 


^1450(1+2) 


-22.45 


fx (0.5-8 keV) 


3.5 X 10"-*^^ ergs 


fx (0.5-2 keV) 


1.8 X 10"-^^ ergs s~^ 


Lx (0.5-2 keV) 


4.5 X 10"^"^ ergs s 



^calculated using HEASARC tool nh 

^AB magnitudes in a 1.4 x FWHM aperture 

■^Rest Frame / band AB magnitude, assuming the 
SDSS average quasar spectrum 



